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5 ' Abstract 

qq ■ Taking into account the characteristics of a free scalar field in elliptic coordinates, a new dynam- 






ical variable is found for the free electromagnetic field. The conservation law associated to this 
variable cannot be obtained by direct application of standard Noether theorem since the symmetry 
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generator is of second order. Consequences on the expected mechanical behavior of an atomic 

^ . system interacting with electromagnetic waves exhibiting such a symmetry are also discussed. 
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I. INTRODUCTION 

Symmetries in physical systems manifest as dynamical constants. In this way, via Noether 
theorem, the homogeneity and isotropy of free space is directly related to the conservation 
of linear momentum P and angular momentum J respectively, and time homogeneity is 
related to energy £ conservation. For a system of particles interacting through intermediate 
fields, global conservation laws are usually translated, via the equations of motion, into local 
conservation laws describing the interchange of dynamical variables between particles and 
fields. Under proper circumstances, the effects of macroscopic materials on the fields are 
synthesized into boundary conditions for the fields. These boundary conditions define a set 
of field configurations (modes) that satisfy them. The mean value of the different dynamical 
variables can be evaluated for each mode inside the region where boundary conditions are 
imposed. If this value of the dynamical variable can change in time just through its flow 
on the boundaries, a conservation law for the dynamical variable associated to the field 
holds. For instance, the geometry of wave guides define the electromagnetic (EM) modes 
and the parameters which characterize them are linked to dynamical properties of the field. 
For rectangular symmetry, the cartesian wave vector k is related to the momentum P of 
the EM field, the frequency to the energy £ per photon and the polarization to the spin 
angular momentum along k, S z . For cylindrical symmetry the energy £ per photon, linear 
momentum P z , orbital momentum along the symmetry axis L z and helicity S z , define the 
parameters that characterize the electromagnetic Bessel modes: frequency u, wave vector 
component along the symmetry axis k z , azimuthal integer m and polarization a . As shown 
twenty years ago, Bessel modes [l|, as well as other electromagnetic configurations coinciding 
with modes inside wave guides of a given geometry, can be generated approximately in free 
space by interferometric means. This has lead to the possibility of creating, in the quantum 
realm, photons with a variety of quantum numbers which, in fact, can be entangled [2|. 

The purpose of this work is to study the dynamical properties of electromagnetic waves 
in elliptic-cylindrical coordinates. The corresponding modes are known as Mathieu fields. 
The cylindrical symmetry leads, under ideal conditions, to propagation invariance along the 
symmetry axis. Four of the five dynamical variables behind the parameters that characterize 
Mathieu modes can be directly identified. The symmetry generator behind the fifth param- 
eter is trivially obtained from the wave equations but its relation to the electromagnetic 



dynamical variable is not, since the symmetry generator is of second order. We shall give 
an explicit expression for this dynamical variable. 

The zeroth order Mathieu beams were first generated in free space by an annular slit 
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illuminated with a strip pattern produced by a cylindrical lens |3j; higher order Mathieu 
beams have already been generated by holographic means [4|. Given the increasing use 
of light to control the motion of atomic systems and microparticles we shall also make 
an analysis of the behavior of an atom in a Mathieu field. Emphasis will be given on 
the dependence of mechanical effects on the values of that parameter and the associated 
dynamical variable. 

II. A MASSLESS SCALAR FIELD IN ELLIPTIC COORDINATES. 



Elliptic-cylindrical coordinates are defined by the transformations 



H 



x + iy = f cosh(u + iv), w G [0, oo), vE[0,2ir) 

z = z, (1) 

where the real valued constant / is half the interfocal distance and the coordinates u and 
v are the so called radial- and angular-like variables, while z is the axial variable. Unitary 
vectors related to a given coordinate x will be written e x . The relevant scaling factor for 
this coordinate system is 

h = h u = h v = /y (cosh2u — cos2w) /2. (2) 

Considering the solution of the wave equation under the assumption of propagation invari- 
ance along the z-direction, 

V 2 ^ = <9 C 2 ^, V(r,t) =iP(f ± )e t ^ z - UJt \ (3) 

Helmholtz equation is obtained — partial derivatives with respect to a given variable x are 
compactly denoted by d x , also the notation d ct = -dt will be used. In elliptic coordinates, 
Helmholtz equation takes the form 

^(aS + fi» + #}iKu,iO = o, kl = ^-e z . (4) 



This equation is separable, tp(u,v) = U(u)V(v), and yields the set of differential equations 

{ajj-6 + 2gcosh2u}J7(u) = 0, (5) 

{d 2 + b-2qcos2v}V(v) = 0, (6) 

known as the modified and ordinary Mathieu equations, in that order, which will be called 
radial and angular Mathieu differential equations from now on. The real constant q is related 
to both half the interfocal distance / and the magnitude of the perpendicular component of 
the wave vector k±_ by 

q = (fk ± /2f . (7) 

From a field theoretical point of view, q is directly related to the perpendicular momentum 
carried by the wave ip. For a given value of q, the possible values of b compatible with the 
boundary condition V(v + it) = V(v) or V(v + 2tt) = V(v) are called the characteristic 
values p, lZJ] - They are usually ordered in ascending values and renamed a n (b n ) for even, 
p = e, (odd, p = o) solutions. The parity of the order n determines if the function is ir- 
on: 27r-periodic for even or odd order n respectively. The mathematical set {ijj(u,vY p ' n ' q ^} is 
complete and orthogonal [8(. 

From Mathieu radial (J3j) and angular ([6]) equations, it is straightforward to construct an 
operator that shares eigenfunctions with the squared transverse momentum, 

B if)(u,v) = b i/;(u,v), 
P 



{ cos 2v dl + cosh 2u d 2 v ) . (8) 
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A physical interpretation of the dimensionless eigenvalue b can be found writing the operator 
B in cartesian coordinates 

x 2 - — ) dl - \ y 2 + — ) d 2 x + 2xyd x d y + xd x + yd y 



^=\ {Iz+k- + h-h+} ~ y Vi = l [z+ l z - } - ^-Vi, (9) 

where the operator l z ± = — % [(r ± fe x ) x V] z , in the context of the quantum mechanics of 
a particle, is proportional to the z— component of the angular momentum with respect to 
either focii of the elliptic-cylindrical coordinate system. 



Previous studies of Helmholtz equation, Eq. (J4]), already report this identification |9l.ll0j|. 
That is the case of the analysis of the spectra of a quantum elliptic billiard for which B 
commutes with the free particle Hamiltonian [9|, LL3J- In fact, j3 = l z +l z - + 2mf 2 H is a 
constant of motion of the classical analogue of this system with m the mass of the particle 
and H the classical Hamiltonian within the billiard. It has been found [9|, [l2| that, if 
l z +l z - > then j3 > 2mf 2 E and the classical trajectory of the particle repeatedly touches 
an ellipse characterized by coshuum = (5 jimp E. While, if l z +l z - < the trajectory always 
lies between the focii touching an hyperbola determined by cosiWi = f3/2mf 2 E. In that 
case f3 is a positive number and l z +l z - has a lower limit given by —2mf 2 E. 

The normalized amplitude of the scalar wave function ip(u,v) is illustrated in Fig. 1 for 
positive and negative values of l\ z+ l z _]. Notice that, in analogy to the classical particle prob- 
lem, for l[ z+ l z _] > ellipses can be observed where the amplitude is approximately constant 
while for l[ z+ l z ~] < hyperbolic patters of similar amplitude are clearly distinguished. 
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Figure 1: (Color Online) Normalized even Mathieu's functions a) ip(u,vp p=e ' n=9 ' q=10 > b) 
ip(u, v)( p=e ' n=9,q=80 > . The corresponding values of the parameter b are 81.6283 and 124.1067, so 
that l[ z+ l z ~] = 61.6283 and h z+ l z _] = —44.1067 respectively. The half focal distance establishes 
length units, / = 1. 



III. DYNAMICAL VARIABLES OF AN ELLIPTIC ELECTROMAGNETIC 
MODE. 

Given a complete set of scalar solutions {^/ K } of the wave equation, a joint set of complete 
electromagnetic modes can be obtained by considering the scalar functions \1/ K as Hertz 
potentials [13|, ll4j. In Coulomb gauge, any solution of the wave equation for the vector 
electromagnetic field A(r, t) 

V 2 A{f,t) = d 2 ct A{f,t) (10) 

can be written as a superposition of modes 

A K (f, t) = A^d ct V x (e 2 vD K ) + _4tp0 [v ± (v • e z ^ K ) - e z V 2 ± ^ K ] (11) 

where k denotes the labels that characterize a given scalar solution \1/ K . The constants A\ 
and A K are proportional to the amplitude of the transverse electric (TE) and transverse 
magnetic (TM) EM fields as can be directly seen from their connection with the associated 
electric E and magnetic B fields, 

E = -d ct A, B = V x A. (12) 

As usual, although expressions may involve complex functions, just the real part of them 
define the corresponding physical quantities. In the case of propagation invariant electro- 
magnetic fields in elliptic coordinates 

A K (f, t) = A (TE) M K + ^(™)N$ K (13) 

where the vector operators are given by the expressions 

M = ^d ct (e u d v - e v d u ) , N = -d z (e u d u + e v d v ) - e z V\. (14) 



As a consequence, 



E K (f,t) = -A {TE) d ct MM K - A {TM) d ct m K , 

B K (r, t) = A iTE) d ct m K - A<- TM) djM K . (15) 



In general, it is expected that the electromagnetic modes given by Eq. (flTl) inherit symme- 
tries of the scalar field with analogous dynamical variables. For the elliptic-cylindrical case, 
invariance under space reflection with respect to the Y axis leads to parity conservation. 
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Meanwhile, invariance under spatial translation along the main direction of propagation of 
the mode is reflected in the fact that the field momentum-like integral 

PJ 4 ' K ' K,) = -*- J d 3 x \(E^ x B®) z ] i = TE, TM (16) 

integrated over the whole space is independent of time. Time homogeneity implies that the 
energy-like integral 

£ {W ' ] = ^j d ' x [^ tf • ^ (?> *) + B { : ] (r, t) ■ B® tt tj\ (17) 

is also constant. In fact, p^ ,K,K and £(*> K > K ') are proportional with ck z /uj the constant 
of proportionality. By construction hk± = K\Juj 2 /c 2 — k 2 would yield the magnitude of the 
transverse component of the momentum which determine the separation constant q, Eq. ((7j). 
Standard quantization rules require a proper normalization selection for the EM modes 
so that each photon carries an energy frw. Using Eq. (12a) from Ref. 151 ] 

/oo /*oo rlir 

dz du dvh u h v ^> K (u,v,z)^' K (u,v,z) = 2'K 2 f 2 s K 5(k z -k' z )5(q-q l )5 ntn , (18) 
■oo JO JO 



_ Vs,2n, q (0)Ve,2n, q (ft/2) Ve,2n+l,q (O^n + l.g W 2 ) 

b e,2n,q ' .(2n) > *e,2n+l,g " 1/2 ,(2n+i) 5 

_ Vo,2n+2, q (0)V:,2 n+ 2, q (*m „ ^l^W^/g ^^ 

6o,2n+2,g — B (2n+2) ) *o,2n+l,g " " i/2 B (2»+i) ' 

where Am and Bm are the standard Mathieu coefficients |6j, |7J , this normalization is trivially 
performed. Defining now the generalized number operator. 

W = \ (aJP+aJP + a«a»t) , [a«, ajf] = 5 id 8^ (20) 

the quantum energy and the momentum along z operators take the form: 

S = J2^^\ V Z = Y,M Z NV, (21) 

allowing the identification of hk z and Hlj with the photon momentum along z and the photon 
energy, in that order. 

For scalar fields and in the case of spacetime continuous symmetries, the generators 
of infinitesimal transformations become good realizations of the corresponding dynamical 
operator. In that sense, the rotation like operator B can be related to the product of 
angular momenta l[ z+ l z -] and the eigenvalue equation B ^ u ,k z ,p,n = b ^ui,k z ,p,n is interpreted 



as a manifestation of the scalar wave function i&w,k z &,n carrying a well defined value of that 
angular momentum product. 

The standard procedure to find the electromagnetic analogue of /[ 2+ / 2 _] would be to apply 
Noether theorem to the field Lagrangian 

C EM = (l/4)(«Vt - d v Ap){d>>A v - d v A») (22) 

for the transformation generated by B on space variables and on the electromagnetic fields 



A^. Standard Noether theorem 16J concerns first order differential operators as generators 
of continuous symmetries, so that, if under an infinitesimal transformation that modify the 
coordinates and field functions by 

Sx, = J2 Xtfu v , 5 A, = J2 <P"M, (23) 

V V 

the Lagrangian is left invariant, the current 

3C 

p = " E ai-to* ~ E a ^ x p) ~ cx p ( 24 ) 

A a 

has zero divergence d v Q v p = 0. As a consequence Q p can be considered as the density of a 
dynamical variable whose integrated value over a volume can change only due to the flux of 
the current 0* through a surface. For the circular cylindrical problem, the assumption of 
isotropy of space via the effects of the infinitesimal rotation generator [f x V] 2 on A^ and 
Xy leads to the identification of 

J z = 2r c [ E^[ fx ^ A ^ 3x +4^/ (ExA) z d 3 x, (25) 

as the z-component of the angular momentum of the electromagnetic field in a volume V 
[lq |. The first integral involves the antihermitian differential operator [r*x V] 2 = d v , and 
is associated to the orbital angular momentum (notice that the standard hermitian angular 
momentum operator is L z = —ihd v ). The second integral is independent of the choice of 
origin, arises from the field variation 5Ay, and is directly related to the polarization of the 



field. It has been identified with the field helicity 17|, ll8|. It can be shown that in the 



Coulomb gauge [l9| 



J= fx(Ex B)d 3 x <b E\(fx A\ ■ ds 

4tvcJ v 47rcJ s lv j 



where S is the surface enclosing the volume V. 

Since the electromagnetic field A^ has a well defined transformation rule under rotations 
which is independent of the origin of space coordinates , the second term in Eq. (1251) . 

S z = -i- / (E x A) z d 3 x, (26) 



is the fourth dynamical variable associated to the Mathieu electromagnetic field. This can 
be directly verified by substituting the general expression for the vectors E and A in terms of 
the elliptic modes. As usual, for a given value of mode indices k, the helicity S z is different 
from zero only if the amplitudes A^™^ and A^ TE ^ are complex. The concept of circular 
and linear polarization of Mathieu waves is analyzed in Ref. |2(| classically. The quantum 
analysis can be carried out in complete analogy with the study in Ref. [18J for Bessel fields 
so that 

s* = h E tr <$> B)ta - M) ~ a - B)a - M)t ) • w 

K 

In the problem treated here, the generator of the transformation B is an hermitian second 
order differential operator obviously related to the isotropy of space. It depends on the 
position of the two focii of the elliptic transversal coordinates and in that sense should be 
analogous to L z . The proposal is to identify the electromagnetic dynamical variable related 
to B with the integral 

B = / J^EitBAi d 3 x, (28) 

41TCJ V ^ 

i 

in complete analogy with Eq. (j25l) . 

In order to corroborate that B is the fifth dynamical variable directly associated to Math- 
ieu electromagnetic waves, notice that 



J2 A?MA* = b A K ' ■ A K 



+ 



v-c-^vi^ 



(29) 



k'kA^A^ + Kk z A^ M) Ar\ 

where the vector C = -1 (M\l/ K /) (fx N_l^ k J + ^ k /N_l^ k — with Nj_ = \d z (e u d u + e v d v ). 

Evaluating the integral in Eq. (1281) over a volume V, using Eq. ( 1291) . the first resulting 
term is proportional to bE K ■ A K involving the same integrals appearing in the energy-like 
density, Eq. ( IT7j) . The second term defines a flux of B through the surface around the 
integration volume V. The third term adds up to the last term in the expression of operator 
B, Eq. ©. 



Eq. (1291) supports the identification of B as the electromagnetic dynamical variable linked 
to the generator B. 

In the quantum realm the corresponding operator is written 

B = Y, n (b + 1) N^ + h(b+ Spj Ni™\ (30) 

In the paraxial limit, the helicity-like factor k z c/cu ~ 1 so that B ~ ^ K h(b + l)NJt . 

Notice that the dynamical variable B for a photon has units of H, although for material 
particles the quantum variable associated to l z +l z - has as natural unit h 2 . The interpretation 
of a dynamical variable for the EM field is usually linked to the interchange of this mechanical 
variable with charged particles or atomic systems. The measurement of B is expected to be 
related to changes in values of l z +l z - although they have different units. It is thus essential 
to clarify how the absorption and/or emission of a Mathieu photon by a particle alters its 
motion. 

IV. MECHANICAL EFFECTS ON ATOMS. 

Theoretical and experimental analysis on the interaction between light and microscopic 
particles, as well as analysis on the interaction of light and cold atoms, have yield very 
important results in the last three decades. In these areas, the use of structured light 



beams with peculiar dynamical properties plays an increasingly important role 2l|. In the 



particular case of Mathieu-like beams, it is possible to generate an elliptical orbital motion 



of trapped microscopic particles [22[ . A detailed theoretical description of this phenomenon 
requires the exhibition of a clear link between the observed motion and the parameters that 
characterize the beam, which are directly related to the mechanical properties of the field 
described here. 

The mechanical effects of a Mathieu electromagnetic wave on a cold atom shall briefly be 
described under the assumption that the atom kinetic energy is low enough to be sensitive 
to light forces but large enough to admit a description in terms of Newton equations. The 



standard semiclassica 



approach is taken, as in the pioneering works by Letokhov 23J and 



Gordon and Ashkin 24J. In this approximation, a monochromatic electromagnetic wave 
describable by a coherent state couples to the dipole moment of an atom. This dipole 
moment /2 12 is related to the electromagnetic transitions between the atom levels that, for 
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simplicity sake, will be taken to have just two accessible options. The gradient of the coupling 
g = ipL\2-E /h = \g\e l( ^ determines the explicit expression for the average semiclassical velocity 
dependent force through the vectors a = V log(|g|) and (3 = V0 24J. The nonlinear Newton 
equations for Mathieu waves have a rich structure that deserves a deep study on its own. 
Here, just some results that illustrate the relevance of the parameter b are reported. Atomic 
transitions with changes in the atomic internal angular momentum 8m = ±1 are proportional 
to e± • E(f,t). For Mathieu waves, this factor is proportional to (ou/c) A TE ± ik z A™ 
reinforcing the interpretation of the latter expression as a signature of circular polarization. 
In order to induce transitions with <5m = ±1 with equal probability, and to avoid any 
complication arising from a component of the electric field along z, it will be assumed that 
A™ = 0. Under these conditions, 

-d v ^d 2 uv ^ + d u ^d 2 j! sinh 1u 



a 






Cy 



(d u ^) 2 + (d v ^) 2 cosh 2u - cos 2v 
h I {d u ty) 2 + (d v V) 2 cosh 2u- cos 2v 



(31) 







d^dl^ - d v ^d 2 J> 



sin2t> 



h I (d u ^) 2 + (d v ^) 2 cosh 2u - cos 2v 






d v vd 2 ,y - d u m 2 ^ 



sinh 2u 



(32) 



(d u ^) 2 + (d v ^) 2 cosh2w-cos2u- 

The expression for the average semiclassical velocity dependent force [24J valid for both 
propagating and standing beams is: 



hT 



(f) 



[{v-a){l-p){l+p)- l + r/2]p+[(v-p)-5uj]a 



(33) 



(1 - p')p'~ l T + 2v ■ a[l -p/p' -p](l +p)~ l 

with T the Einstein coefficient, V = 4k 3 \fli2\ 2 /3h, 5u the detuning between the wave fre- 
quency oj and the transition frequency u; , 5u = lu — uj , p = 2\g\/ ({T /2) 2 + 5u 2 ) a parameter 
linked to the difference D between the populations of the atom two levels, D = 1/(1 + p), 
and finally p' = 2\g\ 2 /\~f'\ 2 , with 7' = (v- a)(l - p)(l + p)' 1 = T/2 + i(-6w + (v- (3)). 

This brief study is focused on the red detuned far off resonance light case so that noncon- 
servative terms arising from the velocity dependence of the force are not dominant 
regime is particularly relevant in the context of optical lattices 26(. Following Ref. 
laser beam is considered with a 67 nm detuning to the red of the 5 2 S\/2 - 5 2 P\/2 transition 
at 795 nm; a 6-10 5 W/cm 2 irradiance is assumed. The trajectories of the atoms are described 
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taking the laser wavelength as unit of length and as unit of time the inverse of the Einstein 
coefficient T which is 3.7 xl0 7 s _1 for the state 5 2 P 1 / 2 of 85 Rb. 

Since Newton equations are highly nonlinear, it is expected that the motion of the par- 
ticles will not have a simple structure. Nevertheless, there are some general characteristics 
of the atom motion that can be predicted without making a numerical analysis. Thus, 
since the z-dependence of the beam is just on its phase, k z z, no longitudinal confinement 
is expected. In the standard paraxial regime, k± << k z , a correspondingly small transfer 
of radial momentum to the atom will occur. Under these conditions, the atom would be 
confined by transversal light potential wells whenever its height is larger than the initial 
atom kinetic energy. Conversely, if k± ~ k z the atom may dwell around different transverse 
wells. Typical trajectories are shown in Figs. 2-3 31]. 
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Figure 2: (Color online) Trajectory of an atom driven by an even Mathieu beam of order n = 
in the paraxial regime with half focal distance / = A, and beam parameters q = 0.0984 and 
b = —0.0048. The initial conditions for the atom are uq = 1.5A, vq = vr/4, zq = 0A, uq = 0.1AT, 

v = -o.ooir, zq = o.ooiAr. 



Using several numerical simulations 28( , the correlations between the time average value 
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Figure 3: (Color online) Trajectory of an atom driven by an odd Mathieu beam or order n = 1 with 
half focal distance / = A, and beam parameters q = 3.5531 and b = —3.5924. The initial conditions 
for the atom are Mo = 1.5A, vq = n/4, zq = OA, uq = O.lAr, vq = — 0.001T, zq = O.OOlAr. Notice 
that the particle dwells on bright zones of the field. 



of the atomic 



LL 



(l z+ l z -)t = m 2 ([(f+ fe x ) x v\ z [(f- fe x ) x v\ z ) t 



and the b value of light mode were studied. 

As illustrated in Fig. 4, in general, the variable l z +lz- exhibits large fluctuations in an 
initial transitory stage. For paraxial beams, there is a time T such that the time average 
l_i_l_ over any interval to < t < to + T becomes independent of to whenever to > T. The order 
of magnitude of T is typically 10 4 r. A rich structure, consistent with the frequent atomic 
recoils due to the light potential wells, could be observed at lower scales. As expected, the 
specific numerical results depend on all the involved variables. Once q and / are fixed, and 
for the same atomic initial conditions, a monotonic nonlinear increase of 1 + 1_ as a function of 

as illustrated in Fig. 5. For non paraxial modes the time average 



b is observed in general ! 



1 + 1_ is highly dependent on the time interval considered due to the persistent increase of 
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the radii of the atom motion illustrated in Fig. 3. 

.JO 
x 10 




x10 



Figure 4: (Color Online) Typical behavior of l z +lz- as a function of time (black)and its average 
(white) for paraxial beams. The initial conditions for the atom are: uq = 1.2A, vq = vr/4, zq = 0A, 
iiQ = 0.2Ar, vq = — O.OOir and zq = O.OOlAr. The odd Mathieu beam is of order n = 7 with half 
focal distance / = A and parameter q = 0.1. The inset graphic is a zoom focused on the initial 
behavior of l z+ l z -. 

V. CONCLUSIONS 



The electromagnetic modes with elliptic-cylindrical symmetry are characterized by their 
polarization and the extended set of parameters k = {u, k z ,p, q, b}, that is the field frequency 
uj and wave vector axial component k z related to the energy £ and z component of the 
linear momentum P z ; the parity p of Mathieu functions and the parameter q which, as in 
the scalar case, is related to the perpendicular component of the wave vector in units of the 
focal distance. It has been proposed and shown that the transformation generator B arising 
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Figure 5: (Color Online) Average 1 + 1_ of the atom, in units of A 2 T _1 , as a function of the b 
parameter of an odd Mathieu beam of order n = 1,3,5,7 with half focal distance / = A and 
parameter q = 0.1. The initial conditions for the atom that yield 1+1_ within the range plotted 
here are u G [.9, 1.5]A, v € [vr/16, vr/4], u € [.1, .2] AT, i) = -O.OOir, i = O.OOlAr. 

from elliptic symmetry is related to the dynamical EM variable B of the electromagnetic 
field, giving a physical significance to b. 

It has been exhibited that the motion of cold atoms in Mathieu beams can be used to 
"measure" B since a strong correlation between the particle product of angular momenta 
1 + 1_ and the parameter b can be approximately isolated by using paraxial TE modes in the 
far off resonance regime. 

The use of light beams with elliptical-cylindrical symmetry has potential applications 
in controlling the mechanical motion of atomic systems which could include nanoparticles. 
Nowadays, it is well recognized that laser-driven nanoparticles have a variety of uses in 
nanofluidics, nanobiotechnology, and biomedicine. However, although the possibility of op- 
tically trapping gold nanoparticles was demonstrated in 1994 30(, developing tweezers for 
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nanoparticles is not straightforward. The gradient forces with conventional beams, fall off 
with particle size. The geometry of elliptical beams adds the focal distance 2/ as parameter 
to tailor gradient forces besides opening the possibility of using the mechanical variable B 
to select over a wider kind of motions. 

The inclusion and further study of this new dynamical variable could also help elucidating 
modern concerns on the "twisted" properties of light, such as why Mathieu functions are 



intelligent states for the conjugate pair exponential of the angle-angular momentum 271 ]. 
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